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Abstract: The application of dynamic modelling for designing surface coal blasting patterns is presented. By 
combining the terms of open cast blasting design pattern, production planning, noise, vibration, fly rock, and 
estimation of drilling and blasting costs, as a decision making tool for planning surface coal blasting 
efficiency, reducing the environmental impacts of noise, vibrations, and fly rock from the surface coal blasting 
can be created. 

The results of the model can be used to design surface coal blasting patterns following mine 
production planning and are based on the theory of the open cast blasting as well. The developed model is 
incorporated in the assessment of drilling and blasting costs, it also makes it easy to compare the results of the 
design in terms of the underlying model for the design of the surface coal blasting. Therefore, this model is 
one of the alternatives that can be used to support the decision making in planning of surface coal blasting 
design patterns effectively. 
Key Words: Surface coal Blasting Design / System Dynamics Model / Cast Blasting Model 



1. INTRODUCTION 

Drilling and blasting for overburden and coal in open cast mines is estimated for 55% of total mining 
costs in open cast mining, for this reason the planning and design of drilling and blasting to minimize this cost 
are necessary [1]. The selection of optimum burden, or powder factor, are an importance parameters in open 
cast blasting design pattern. The interest in this type of blasting is based on a reduced investment in 
machinery, operation and maintenance, as there is less manipulation of material due to the fact that a large 
volume of rock, from 40 to 60%, can be projected [2]. 

Not only blasting design parameters to be decided but environmental impact parameters from blasting 
also be considered. Thus, the tool that can include both groups of parameter is important. 

A dynamics model of surface coal blasting design (DMSCBD) is followed the theory of open cast 
blasting and is included: the equations and design criteria for surface coal blasting, secondly, environmental 
control equations for controlling the impacts of noise, vibrations, and fly rock in blasting [3-4], and finally, 
cost estimation from the drilling and blasting design. This model is created by Vensim Software, which is 
generally used in system dynamics modelling design [5]. The aim of DMSCBD is to help the designing of 
optimum patterns, safely charge of explosive per delay, and optimum cost per production. 

2. SURFACE COAL BLASTING DESIGN 

Cast blasting is being applied with great success in coal mines in South Africa, United States, Australia 
and Canada [2]. The method may be defined as the use of explosive for the purpose of fragmenting and 
providing displacement of the overburden to the final spoil pile. The maximum forward displacement is the 
primary requirement in this blasting method [3]. 

In case of surface coal blasting design, while using open cast blasting, some methodology or equation 
was presented. In 1980, D' Appolonia Consulting Engineers of Pittsburgh developed a set of nomographs for 
designing blasting rounds and mines. Blasting data from actual blasting practices in approximately 100 surface 
coal mining operations were collected to develop these nomographs. The nomographs are shown in Figs. 1, 2, 
3 and 4 [1-2]. Paul et al. (1987) [1] developed a computer aided blasting-round design in an arctic coal mine 
using the nomographs. It was developed in Fortran 77. 

The nomographs contain 3 sets of variables: input variables, intermediate variables and output variables. 
Input variables are: the diameter of the blast hole (D e ), strain energy factor (£), bench height (H), desired 
distance the blasted material is to be thrown (R), and density of explosive (P). The intermediate variables are 



those which are calculated in the process of determining output variables and they consist of 5 constants as Cj, 
C 2 , C 3 , K] and K 2 , blastability factor (E ), hole loading length {Hj), hole loading density (P cl ) and powder 
factor (Q). The output variables are those which are used for production drilling and blasting design. They 
include the burden of the blast hole pattern (B, Bj), spacing between the holes {S), stemming length (S t ), and 
total charge weight per hole (W). The following empirical equations are the basis of the nomographs: 
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The optimum burden B opt can be found while B=Bj by varying K] and K 2 . D' Appolonia uses a rule of 
thumb which is K 2 = Kj to find B opt [2]. 

The relationship between compressive strength of rock (CS) and strain energy factor (E) is shown in 
Table 1 [2]. 
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Other equations for the calculated patterns of surface coal blasting followed a prototype dynamics model 
of bench blasting design [6]. 
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Fig 1 Nomograph 1 



Fig 2 Nomograph 2 
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Fig 3 Nomograph 3 



Fig 4 Nomograph 4 



3. SYSTEM DYNAMICS 

System dynamics is an approach to understanding the behavior of complex systems over time. It is a 
powerful methodology and computer simulation modeling technique for understanding, and discussing 
complex issues and problems [7]. 

System dynamics which founded by Prof. J.W. Forrester in 1950 [8], is a theory of system structure and a 
set of tools for representing the structure of complex systems and analyzing their dynamic behavior. In Fig 5 
shows the generic structures for creating a dynamics model. 
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Fig 5 Generic structure of system dynamics model 

P. Sontamino and C. Drebenstedt (2012) [6] developed a prototype dynamics model of bench blasting 
design using Vensim software. The prototype model covered bench blasting design parameters, cost of drilling 
and blasting parameters and environmental impact parameters. It was an alternative and flexible tool to 
support user to design bench blasting patterns and decided the suitable condition by terms of economic and 
environmental control. 



4. METHODOLOGY 

In this paper, model structures and equations are created in Vensim software [5] following system 
dynamics theory and open cast blasting design pattern by using the equations converted from the nomographs 
[1-3]. The model is developed based on the prototype dynamics model of bench blasting design [6]. By 
changing some parameters and equations from the prototype dynamics model of bench blasting design, the 
dynamics model of coal blasting design pattern can be determined. The structure of the model is shown in Fig 
7. 



The simple flowchart of DMSCBD is presented in Fig 6. 
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Fig. 6 Simple flow chart of DMSCBD 



5. LIST OF VARIABLES 



There are 30 main input variables, 34 main output variables, and 4 main criterion variables in Table 4, 5 
and 6 respectively. However, in the structure of the model, it is necessary to add more variables than are 
shown in the table below, such as function for selecting values of coefficient factors in the equations, the 
converter variables for correcting unit of variables, etc. 

Every variable is defined as a symbol in a short description which will be used in the model structure. 



Table 4 List of main input variables 


in DMSCBD 


Name 


Symbol 


Value 


Ref. 


1 . A value of ground 
vibration equation 


A 


0.607 

(Sandstone) 


[4] 


2. Alpha value of ground 
vibration equation 


Alpha 


1.47x10 5 
(Sandstone) 


[4] 


3 . B value of ground 
vibration equation 


B 


1.463 

(Sandstone) 


[4] 


4. Blast hole inclination 
(degree) 


BHI 




(Vertical) 


* 


5. Blasting production 
planning (tons/day) 


BPP 


3,200- 
27,200 


Sc 


6. Cap per hole (cap/hole) 


CH 


1 


=!= 


7. Compressive strength 
(Mpa) 


CS 


1,5, 10 


Sc 


8. Hole Diameter (inch) 


D 


5,7,9 


Sc 


9. Density of ANFO 
(kg/cu.m.) 


DANFO 


800 

(Loose) 


[9] 


10. Density of Fuel Oil 
(kg/Liter) 


DFO 


0.8 


[3] 


11. Density of High 
Explosive (kg/cu.m.) 


DHE 


1,200 
(Emulsion) 


[3] 


12. Density of materials 
(kg/cu.m.) 


DM 


1,600 

(Soil) 


[3] 


13. Input drilling rate 
(m/min) 


ID 


0.35 


[3] 


14. Input bench high (m) 


H 


9, 11, 13 


Sc 


15. Bottom charge input (m) 


ILf 


1 


* 


16. Subdrilling input (m) 


J 





* 


17. K value of ground 
vibration equation 


K 


713.23 
(Sandstone) 


[4] 



Table 5 List of main output variables in DMSCBD 



Name 


Symbol 


Name 


Symbol 


1 . Cumulative 
blasting 
production (tons) 


CBPP 


2. Optimum 
burden (m) 


B-opt 


3. Optimum 
spacing (m) 


S-opt 


4. Optimum 
stemming (m) 


St-opt 



Name 


Symbol 


Value 


Ref. 


18. Coefficient of maximum 
distance of fly rock 


Kfr 


260 


[3] 


19. Coefficient of maximum 
diameter of fly rock 


Kfs 


0.1 


[3] 


20. Reference Pressure (bar) 


P0 


2.0x10"'° 


[3] 


21. Peak particle velocity 
(mm/s) 


PPV 


8 

(DIN 4150) 


[3] 


22. Distance between 
blasting area and measuring 
area (m) 


R 


1,500 


* 


23. Distance to be thrown 
(m) 


Rt 


1,5, 10 


Sc 


24. Ratio of AN (%) 


RAN 


94.5 


[9] 


25. Unit cost per kg of AN 
(euro/kg) 


UAN 


0.42 


q 


26. Unit cost of Cap 
(euro/cap) 


UC 


0.62 

(Electric) 


q 


27. Unit cost of drilling 
(euro/hr) 


UD 


35 


* 


28. Unit cost per Liter of 
Fuel Oil (euro/Liter) 


UFO 


0.67 


q 


29. Cost of High Explosive 
(euro/kg) 


UHE 


2.29 


q 


30. Work time for drilling 
(hr/day) 


WTD 


8 


* 



* = Assumption value; 

Sc = variables use of simulation in 

q = in the query 



scenarios. 



Name 


Symbol 


Name 


Symbol 


5. Charge Zone 
(m) 


CZ 


6. Blast hole 
length (m) 


1 


7. Breakage 
volume per hole 
(cu.m.) 


VR 


8. Breakage 
weight per hole 
(tons) 


BWH 



Name 


Symbol 


Name 


Symbol 


9. Weight of 
ANFO (kg) 


Qc 


10. Total Weight 
of Explosive (kg) 


Qb 


11. Weight of 
High Explosive 
(kg) 


Qf 


12. Number of 
holes to drill 
(hole/day) 


NH 


13. Time for 
drilling (hr/day) 


TD 


14. Number of 
drilling machine 
required 


NDM 


15. Cost of High 

Explosive 

(euro/day) 


CHE 


16. Cost of AN 
(euro) 


CAN 


17. Cost of 

ANFO 

(euro/day) 


CANFO 


18. Cost of Fuel 
Oil (euro) 


CFO 


19. Cost of Caps 
(euro/day) 


CC 


20. Total cost of 

explosive 

(euro/day) 


TCE 


21. Total cost of 

blasting 

(euro/day) 


TCB 


22. Total cost of 

drilling 

(euro/day) 


TCD 


23. Total cost of 
drilling and 
blasting 
(euro/day) 


TCDB 


24. Cumulative 
cost of drilling 
and blasting 
(euro) 


CCDB 


25. Cost per 

production 

(euro/ton) 


CPP 


26. Q max (kg) 


Qm 



Name 


Symbol 


Name 


Symbol 


27. Maximum 


MHpD 


28. Maximum 


MDFR 


hole per delay 




distance of Fly 




(hole) 




Rock (m) 




29. Noise from 


db 


30. Estimate air 


P 


Blast (dB) 




pressure from the 
blast (bar) 




3 1 . Fly rock 


FRD 


32. Yield of 


RA 


diameter size 




broken rock 




(m) 




(m 3 /m) 




33. Decision 


FDM 


34. Real 


RP 


making for all 




Production 




blasting design 




(cu.m./day) 




conditions 








(l=yes, 0=no) 









Table 6 List of criterion variables in DMSCBD 



Name 


Symbol 


Recommended 


Ref. 


1 . Powder factor 
(kg/m 3 ) 


PF 


0.3-0.8 


[2] 


2. Vibration problem 
(l=yes, 0=no) 


VP 


MHpD>=l 


=!= 


3. Noise problem 
(l=yes, 0=no) 


NP 


<172dB 


[9] 


4. Fly rock problem 
(l=yes, 0=no) 


FP 


R>MDFR 





: Assumption criteria 



6. MODELS STRUCTURES 

Model structures of DMSCBD are separated into 6 sub-models (Fig 7). 

■ Blasting Design Patterns ■ Drilling and Blasting Conditions 

■ Blasting Production Planning ■ Costs of Drilling and Blasting 

■ Environmental Control ■ Decision Making 

In the model structure, the input variables are shown in O symbol, the output and converter variables are 
shown in normal text, the criterion variables are shown in O symbol, the cumulative variables are shown in 
I I symbol, the rate variables are shown in ^**° symbol, and all of the sub-models are connected by 
variables with < > symbol. 
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[Sub-model: Blasting Design Patterns] 
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[Sub-model: Environmental Control] 
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[Sub-model: Decision Making] 
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Fig. 7 Model Structure of DMSCBD 

7. MODEL SCENARIOES 

From Table 4: List of main input variables in DMSCBD, four variables are taken into account to create 
sample simulation scenarios, including hole diameter (D), compressive strength of rock (CS), Bench high (H), 
and the distance to be thrown (Rt). The aim of these sample simulations is to find optimum burden (B-opt), 
spacing (S-opt) and stemming (St-opt). By changing each variable, there are 81 scenarios to be created 
following Table 7 



Table 7 Metrix variables in the scenario simulations 



Variables 


Valuel 


Value2 


Value3 


1 . Distance to be thrown Rt 
(m) 


1 
(Rtl) 


5 
(Rt5) 


10 
(RtlO) 


2. Bench high H 
(m) 


9 
(H9) 


11 
(Hll) 


13 
(H13) 



Variables 


Valuel 


Value2 


Value3 


3. Hole diameter D 
(inch) 


5 
(D5) 


7 
(D7) 


9 
(D9) 


4. Compressive strength CS 
(MPa) 


1 

(CS1) 


5 
(CS5) 


10 
(CS10) 



In each scenario, the symbol is used to make the short explanation, for instance "RtlH9D5CS10" means: 
the design condition of the distance to be thrown equals 1 meter, bench high equals 9 meters, hole diameter 
equals 5 inches, and compressive strength of the material to be blasted equals 10 Mpa. 

8. SIMULATION RESULTS 

From 81 scenarios, the result of the optimum burden (B-opt), spacing (S-opt), and stemming (St-opt) are 
determined in Table 8. 



Table 



8 Optimum burdens 


, spacings 


No. 


Scenarios 


B-opt 


S-opt 


St-opt 


1 


RHH9D5CS1 


3,72 


4,82 


4,06 


2 


RHH9D5CS5 


3,76 


4,58 


4,02 


3 


RHH9D5CS10 


3,82 


4,30 


3,97 


4 


RHH9D7CS1 


4,22 


7,20 


5,04 


5 


RHH9D7CS5 


4,28 


6,83 


5,00 


6 


RHH9D7CS10 


4,34 


6,41 


4,94 


7 


RHH9D9CS1 


4,56 


9,69 


5,86 


8 


RHH9D9CS5 


4,62 


9,18 


5,81 


9 


RHH9D9CS10 


4,70 


8,60 


5,75 


10 


RHH11D5CS1 


4,06 


4,91 


4,33 


11 


RHH11D5CS5 


4,11 


4,65 


4,28 


12 


RHH11D5CS10 


4,17 


4,36 


4,23 


13 


Rt1H11D7CS1 


4,67 


7,32 


5,43 


14 


RHH11D7CS5 


4,73 


6,94 


5,37 


15 


Rt1H11D7CS10 


4,80 


6,51 


5,31 


16 


Rt1H11D9CS1 


5,12 


9,87 


6,35 


17 


RHH11D9CS5 


5,18 


9,26 


6,29 


18 


RHH11D9CS10 


5,27 


8,68 


6,22 


19 


RHH13D5CS1 


4,36 


4,94 


4,54 


20 


RHH13D5CS5 


4,41 


4,67 


4,49 


21 


RHH13D5CS10 


4,47 


4,38 


4,44 


22 


RHH13D7CS1 


5,05 


7,40 


5,74 


23 


RHH13D7CS5 


5,11 


7,00 


5,68 


24 


RHH13D7CS10 


5,19 


6,57 


5,61 


25 


RHH13D9CS1 


5,58 


9,88 


6,75 


26 


RHH13D9CS5 


5,66 


9,35 


6,69 


27 


RHH13D9CS10 


5,74 


8,77 


6,61 



and 



stemmings in 


31 scenarios 


No. 


Scenarios 


B-opt 


S-opt 


St-opt 


28 


Rt5H9D5CS1 


3,77 


3,93 


3,82 


29 


Rt5H9D5CS5 


3,81 


3,76 


3,79 


30 


Rt5H9D5CS10 


3,86 


3,56 


3,76 


31 


Rt5H9D7CS1 


4,31 


5,86 


4,77 


32 


Rt5H9D7CS5 


4,36 


5,60 


4,74 


33 


Rt5H9D7CS10 


4,42 


5,31 


4,70 


34 


Rt5H9D9CS1 


4,68 


7,85 


5,56 


35 


Rt5H9D9CS5 


4,74 


7,50 


5,52 


36 


Rt5H9D9CS10 


4,81 


7,11 


5,48 


37 


Rt5H11D5CS1 


4,11 


3,97 


4,06 


38 


Rt5H11D5CS5 


4,15 


3,80 


4,03 


39 


Rt5H11D5CS10 


4,21 


3,60 


4,00 


40 


Rt5H11D7CS1 


4,75 


5,95 


5,12 


41 


Rt5H11D7CS5 


4,80 


5,69 


5,08 


42 


Rt5H11D7CS10 


4,87 


5,39 


5,04 


43 


Rt5H11D9CS1 


5,22 


7,94 


6,01 


44 


Rt5H11D9CS5 


5,29 


7,59 


5,97 


45 


Rt5H11D9CS10 


5,37 


7,20 


5,92 


46 


Rt5H13D5CS1 


4,40 


3,99 


4,26 


47 


Rt5H13D5CS5 


4,45 


3,81 


4,22 


48 


Rt5H13D5CS10 


4,50 


3,61 


4,18 


49 


Rt5H13D7CS1 


5,12 


6,00 


5,40 


50 


Rt5H13D7CS5 


5,18 


5,73 


5,36 


51 


Rt5H13D7CS10 


5,25 


5,43 


5,31 


52 


Rt5H13D9CS1 


5,68 


8,02 


6,37 


53 


Rt5H13D9CS5 


5,75 


7,67 


6,33 


54 


Rt5H13D9CS10 


5,83 


7,26 


6,27 



No. 


Scenarios 


B-opt 


S-opt 


St-opt 


55 


RH0H9D5CS1 


3,82 


3,19 


3,60 


56 


RH0H9D5CS5 


3,86 


3,08 


3,58 


57 


RH0H9D5CS10 


3,91 


2.94 


3,55 


58 


RH0H9D7CS1 


4,38 


4,77 


4,51 


59 


RH0H9D7CS5 


4,43 


4,60 


4,49 


60 


RH0H9D7CS10 


4,49 


4,40 


4,46 


61 


RH0H9D9CS1 


4,79 


6,39 


5,27 


62 


RH0H9D9CS5 


4,85 


6,15 


5,25 


63 


RH0H9D9CS10 


4,92 


5.88 


5,22 


64 


Rt10H11D5CS1 


4,15 


3,22 


3,81 


65 


Rt10H11D5CS5 


4,20 


3,10 


3,79 


66 


Rt10H11D5CS10 


4,25 


2,96 


3,77 


67 


Rt10H11D7CS1 


4,82 


4,84 


4,82 


68 


Rt10H11D7CS5 


4,87 


4,66 


4,80 


69 


Rt10H11D7CS10 


4,93 


4,45 


4,77 


70 


Rt10H11D9CS1 


5,32 


6,47 


5,68 


71 


Rt10H11D9CS5 


5,38 


6,23 


5,65 


72 


Rt10H11D9CS10 


5,46 


5,95 


5,62 


73 


RI10H13D5CS1 


4,44 


3,22 


3,99 


74 


RI10H13D5CS5 


4,49 


3,10 


3,97 


75 


RI10H13D5CS10 


4,54 


2,96 


3,94 


76 


RI10H13D7CS1 


5,18 


4,87 


5,08 


77 


RI10H13D7CS5 


5,24 


4,69 


5,05 


78 


RI10H13D7CS10 


5,13 


4,48 


5,02 


79 


RI10H13D9CS1 


5,77 


6,53 


6,01 


80 


RI10H13D9CS5 


5,84 


6,28 


5,98 


81 


RI10H13D9CS10 


5,91 


6,00 


5,94 



From Table 8, the results of optimum burdens, spacings, and stemmings change following 
variables as Rt, H, D, and CS which is shown in summary result (Fig 8). 
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Fig. 8 Summary result of scenarios 

In Fig. 8 shows the relationship of variables such as, when increasing the distance to be thrown (Rt), 
blasting pattern should be reduced. When increasing bench high (H) and/or hole diameter (D), blasting pattern 
increase, finally when the compressive strength of the material (CS) is increasing, optimum burden also 
increase but spacing and stemming are decreasing, etc. 

9. CONCLUSIONS 

A dynamics model of surface coal blasting design shows the results in many scenarios automatically 
when changed input variables. Thus, It can helps to design surface coal blasting patterns following the theory 
of open cast blasting design. 

However, this model is a generic model. It needs to be developed for more user friendly. An adjustment 
and update value of input variables related to the site conditions before used are necessary. 

In the future, DMSCBD can be extended and included variable in another process of coal mining such as 
loading, transportation, crushing and processing which will lead to the more useful tool for a decision making 
in coal mine planning. 
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